Excited hadrons on the lattice: Mesons 
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We present results for masses of excited mesons from quenched calculations using chirally improved 
quarks at pion masses down to 350 MeV. The key features of our analysis are the use of a matrix 
of correlators from various source and sink operators and a basis which includes quark sources with 
different spatial widths, thereby improving overlap with states exhibiting radial excitations. 



PACS numbers: 11.15.Ha 

I. INTRODUCTION 

Ground-state hadron spectroscopy has been at the 
forefront of the studies performed by the lattice QCD 
community. The exponential decay behavior of the mass 
eigenstate contributions to Euclidean-time correlators 
makes extraction of the lightest state straightforward: at 
large time separations, the ground-state mass naturally 
dominates the correlator, which then displays a single 
exponential in time. Excited states, however, appear as 
subleading exponentials, which in practice are rather dif- 
ficult to separate, not only from the ground-state contri- 
bution, but also from those due to other excited states. 

A few methods are currently in use to deal with this 
problem. A straightforward fit to a finite sum of exponen- 
tials is sometimes possible, usually when high statistics 
are available. Constrained fitting [lj and the Maximum 
Entropy Method are also options, but these, like the 
first method, can encounter problems when faced with 
states which lie close in mass, or when unphysical contri- 
butions appear due to quenching (ghosts). 

We employ a different method, one based upon the 
variational method devised by Michael and later 
elaborated by Liischer and Wolff 4]. The distinguish- 
ing characteristic of our approach is the use of several 
hadron interpolators which contain quark wavefunctions 
covariantly smeared to approximate Gaussians of differ- 
ent widths 5]. Using such a basis of source and sink 
operators, it is our plan to improve overlap with physical 
states which may involve a radial excitation. In a recent 
paper [6j , it has been demonstrated that this method also 
clearly separates ghost contributions. 

In the present work, we describe the method in detail 
and report results for the meson sector. Another paper 
will contain our latest results for baryons. 

II. THE METHOD 

Our calculation is based upon the variational method 
0, Q . The central idea is to use several different interpo- 
lators Oi, i = 1, . . . N with the quantum numbers of the 
desired state and to compute all cross correlations 

C(% = (0,(4)^(0)). (1) 



In Hilbert space these correlators have the decomposition 
C(t)ij = ^(0|O l |n)(n|O]|0) e -* A/ " . (2) 

n 

Using the factorization of the amplitudes one can show £| 
that the eigenvalues (4) of the generalized eigenvalue 
problem 

C(4)«W = A (fc) (4)C(4 )^ fc) , (3) 

behave as 

A (fc) (4) oc e~* Mk [ 1 + 0(e -t AMk ) ] , (4) 

where Mk is the mass of the fc-th state and AMk is the 
difference to the mass closest to Mk 0. In Eq. © the 
eigenvalue problem is normalized with respect to a times- 
lice to < t. 

Equation Q shows that the eigenvalues each decay 
with their own mass: The largest eigenvalue decays with 
the mass of the ground state, the second largest eigen- 
value with the mass of the first excited state, etc. Thus, 
the variational method allows one to decompose the sig- 
nal into those for ground and excited states, as well as 
ghost contributions 6], and therefore, simple, stable two- 
parameter fits become possible. 

We remark that also for the regular eigenvalue prob- 
lem, 

C{t)4 k) = A (fc) (4)^ fc) , (5) 

a behavior of the type (J3J can be shown £| . We find that 
in a practical implementation with our sources the results 
and the quality of the data are essentially unchanged 
whether we use the regular or the generalized eigenvalue 
problem. 

The key to a successful application of the variational 
method is the choice of the basis interpolators. They 
should be linearly independent, and as orthogonal as pos- 
sible, while at the same time, able to represent the phys- 
ical state as well as possible. Furthermore, they should 
be numerically cheap to implement. 

For meson spectroscopy it is well known that different 
Dirac structures can be used to construct interpolators 
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TABLE I: Quantum numbers of the interpolators © for 
different choices of V. We remark that the classification with 
respect to C is for flavor degenerate interpolators only. 
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with the quantum numbers one is interested in. In this 
paper we consider interpolators of the form 







(6) 



where fi are flavor labels. We use matrix/vector nota- 
tion for Dirac and color indices and T is an element of 
the Clifford algebra. Depending on what combination of 
gamma- matrices one uses for T, the interpolator O will 
have different quantum numbers. We list these in TableQ] 

However, it is well known (see, e.g., 0, @) that cor- 
relating operators with different Dirac structures alone 
does not provide a sufficient basis to obtain good overlap 
with excited states. For many excited hadrons the spa- 
tial wavefunctions are expected to have nodes. In |5lllCI| 
it was proposed, and tested on small lattices, to use Ja- 
cobi smeared quark sources of different widths to allow 
for nodes in the radial wavefunction. Other lattice ef- 
forts 0, 0| have also seen the need of using spatially 
extended operators. 

In a lattice spectroscopy calculation the hadron corre- 
lators are built from quark propagators D^ 1 acting on a 
source s, 



* * b c * 



(7) 



y,p,c 



If the source is point-like, i.e., s = sq, with 

4 a,< %o), = W,Q)8 pa 6 ca , (8) 

then the two quarks in JfJJl both sit on the same lattice 
site. Certainly this is not a very physical assumption. 

The idea of Jacobi smearing |l3|, [l4| is to create an 
extended source by iteratively applying the hopping part 
of the Wilson term within the timeslice of the source: 

N 



a (a,a) = M fl (a,a) ^ M = ^ K™ H n , 

n=Q 

3 

H(x,y) = [U j (x,0)6(x + j,y) 

+ U 3 {x~hQ)U{x~hv)\ ■ (9) 



Applying the inverse Dirac operator as shown in Q con- 
nects the source at timeslice t = to the lattice points at 




FIG. 1: Profiles P(r) of the narrow and wide source as a func- 
tion of the radius r (for the f3 = 7.90 lattice). The symbols 
are our data points. The curves are the target Gaussian dis- 
tributions (with half width a) which we approximate by the 
profiles P(r). 



timeslice t. There an extended sink may be created by 
again applying the smearing operator M. 

Jacobi smearing has two free parameters, the hopping 
parameter k and the number of smearing steps N. They 
can be used to create sources and sinks with approxi- 
mately Gaussian profiles of different widths. In Fig. 2] we 
show two such profiles P(r) as a function of the radius r. 
For mapping these profiles we use the definition 



P(r) = J2 S 



\y\ 



E 



(£0). 



(10) 



The parameters k and N were adjusted such that the 
narrow source approximates a Gaussian with a half width 
(i.e., standard deviation) of 0.27 fm and the wide source a 
Gaussian with a half width of 0.41 fm. The corresponding 
smearing parameters are listed in Table ^ The values 
of the smearing parameters are chosen such that simple 
linear combinations of our two sources allows one to ap- 
proximate the ground- and first-excited state wavefunc- 
tions of the 3-d harmonic oscillator (with a half width 
of a = 0.33 fm for the ground state). We stress that 
we do not impose such linear combinations, but rather 
use the simple n, w sources for the basis interpolators in 
the variational method and leave it to the simulation to 
determine the physical superpositions. 



TABLE II: Parameters of our simulation. We list the lattice 
size, the inverse coupling (3, the number of configurations, the 
lattice spacing a, the cutoff a -1 , and the smearing parameters 
TV and k for the narrow and wide sources. 
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FIG. 2: Effective mass plots for pseudoscalar (PS) and vector (V) mesons from our coarse (16 3 x 32, a = 0.148 fm, with 
am q = 0.05,0.1,0.2 from top to bottom) and fine (20 3 x 32, a = 0.119 fm, with am q = 0.04,0.08,0.16) lattices. Both ground 
and excited states are shown, along with the M ± a M results (horizontal lines) from correlated fits to the corresponding time 
intervals. The numbers in the plots are the % 2 P er degree of freedom values for the fits to the excited-state plateaus. 



The interpolators which we actually consider in our 
correlation matrix are flavor triplet operators of the form, 

O r = TiTd. (11) 

This avoids the need to calculate disconnected pieces. 
For a given Clifford element T both u and d sources can 
be wide or narrow which gives the 4 possibilities (the 
subscripts n, w denote the type of smearing used) 

^ Km — ^ n ^- ^ n ? Own = ^ w ^ 3 
O nw — U n J^ d w , O ww — U w T*d w . (1^) 

In addition, we have two choices for T in the vector and 
pseudoscalar sectors such that for these cases we have a 



basis of 8 interpolators. For scalars and pseudovectors 
we restrict ourselves to 4 interpolators. 

In the case of degenerate quark masses, 0^ I]n and 0^ w 
give rise to identical correlators. This fact reduces the 
correlation matrices to 6 x 6 for vector and pseudoscalar 
states and to 3 x 3 for scalar and pseudovector states. 
For the strange mesons, obtained by replacing the down 
quark in (|ll(l by a strange quark, the degeneracy is lifted 
and we are free to work with the full 8x8 and 4x4 
correlation matrices. 

Our quenched study is done with the Liischer-Weisz 
gauge action at two different values of the inverse 
gauge coupling (3. The lattice spacing a was determined 
in [lfjj using the Sommer parameter. We use lattices of 
two different sizes, such that the spatial extent in phys- 
ical units is kept constant at 2.4 fm. This allows us to 
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compare the results at two different values of the cutoff, 
a" 1 = 1680 MeV and a" 1 = 1350 MeV. For the fermions 
we use the chirally improved Dirac operator [17| . It is 
an approximation of a solution of the Ginsparg- Wilson 
equation [18j, with good chiral behavior 19]. The pa- 
rameters of our calculation are collected in Table HU We 
determine the strange quark mass via interpolations in 
the heavy-quark mass which match the (light-quark mass 
extrapolated) pseudoscalar K meson mass to the physi- 
cal value. (These configurations have also been used in 
another study [2(J to determine low energy constants.) 

We remark that the CI operator has one term which 
is next-to- nearest neighbor. This has to be kept in mind 
when selecting fit ranges for masses. Exactly what is 
to be considered a safe minimum time separation before 
fitting is not, however, a simple matter. In the present 
case, we limit ourselves to values larger than At = 2. 



III. RESULTS 

A. Effective masses and fit ranges 

Since this is a report on lattice spectroscopy, we be- 
gin the discussion of our results in the standard manner: 
with effective masses. Having performed the necessary 
diagonalizations of the correlator matrices, we determine 
the effective masses from ratios of eigenvalues on adjacent 
timeslices: 



aM 



eff 



= In 



U (fc) (i + 1) 



(13) 



Errors for such quantities are determined via a single- 
elimination jackknife procedure. 

We note that the quality of the plateaus which we ob- 
tain here can depend upon the basis chosen. This is sim- 
ply due to the fact that the various interpolators have 
different overlap with the states which contribute to the 
higher-order corrections in Eq. 10} (e.g., see Ref. [Tl|p. 
An appropriate choice of operator combination (includ- 
ing the possibility of further limiting the basis) can there- 
fore minimize such corrections, improving the plateaus in 
effective mass (and correspondingly, in the eigenvector 
components). 

In Fig. |21 we present our effective mass plots for what 
we consider to be our optimal operator combinations 
(nTn, nj^Tn, wy4Tw, wj^Tw for pseudoscalar and vec- 
tor mesons on our coarse lattices; wjiTn, n^^Tw, wjiTw 
for the same mesons on our fine lattices; and nTn, nTw, 
wTw otherwise). Shown are results for the pseudoscalar 
and vector mesons from three quark masses on both sets 
of lattices. In each case, the effective masses from the 
largest two eigenvalues (largest three for the vectors on 
the fine lattice) are displayed. The horizontal lines in 
the plots mark the M ± a M values which arise from cor- 
related fits to the corresponding time intervals (we fit 
only when we see a plateau of at least three successive 
effective-mass points). Also given in the plots are the 
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FIG. 3: Eigenvector for the first-excited vector meson on our 
20 3 x 32 fine lattices (a » 0.12 fm, am q = 0.06, t = la). The 
eigenvector on each time slice is normalized such that \v\ — 1. 
There is a clear jump in the relative components from t = 3a 
to t — 4a, a plateau forming afterwards. In this situation, we 
fit starting at t = 4a, even though the corresponding effective 
mass plot shows a plateau starting sooner. We apply similar 
restrictions to all our fit ranges. 



X 2 /d.o.f. values for the fits to the excited-state plateaus. 
Here, one can see questionable plateaus, and hence poor 
fits, for the coarse-lattice excited vectors at high quark 
mass, but this is the only place where such high x 2 /d.o.f . 
values occur. For all other fits, x 2 /d.o.f . < 1. 

To be sure that we are dealing (within statistical er- 
rors) with a single mass eigenstate, we check that the cor- 
responding eigenvectors remain constant over the same 
region where we fit the eigenvalues for the mass. An 
example of such a test is shown in Fig. |21 where we dis- 
play the second eigenvector for one of the vector mesons 
on our fine lattices. It is obvious here that plateaus for 
this state cannot be trusted before t = 4a, even though 
when looking at the corresponding effective masses, one 
is tempted to start the fit at t — 3a (see the right-most 
column in Fig. EJ) . We make sure that all our fit ranges 
obey (within errors) such restrictions. (This procedure 
accounts for the enlargement of error bars for the afore- 
mentioned state since we reported preliminary results for 
the mesons [2l[. We are simply being more conservative 
now.) 



B. Quark mass dependence 

In Figs.QHSJ we plot the resulting hadron masses ver- 
sus the ground-state pseudoscalar mass-squared (oc m q ), 
all in units of the Sommer parameter rQ. All plots display 
results from both our coarse and fine lattices. The ver- 
tical lines mark the values of r^M^ which arise when we 
use the strange quark mass in the degenerate-quark-mass 
pseudoscalar; i.e., these mark the point where the mass 
of each of the light (u, d) quarks in the corresponding 
meson equals the strange quark mass. 

Figure 01 shows the excited pseudoscalar masses. The 
results from both lattice sets appear consistent with the 
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FIG. 4: Excited-state pseudoscalar (PS) masses vs. for 
both lattice spacings. The quark masses are degenerate. All 
quantities are in units of the Sommer parameter, ro. The 
diamond represents the experimental point and the vertical 
lines mark the values of r^M^ corresponding to the physical 
strange quark mass. 



and this despite the fact that our method removes the 
dominant influences of ghosts 6] . The entire scalar plot 
also appears to be shifted vertically, making one wonder 
whether the ground-state ud scalar meson even corre- 
sponds to the a (980). Other lattice studies 0,113 find 
that, in fact, it does not. 

Our interpolators may be a poor choice of opera- 
tors for the pseudovector and scalar states. One might 
achieve better results by altering the S'-wave nature of 
our smeared quark wavefunctions. In a future study, 
we plan to apply covariant derivatives to our smeared 
sources, providing interpolators which better mimic P- 
wave orbital excitations and, hopefully, improving over- 
lap with the pseudovector and scalar mesons. 

In Fig. EI we present our results for strange-light pseu- 
doscalar and vector mesons. The abscissa is again TqM^, 
but now this only corresponds to the light-quark mass 
and we therefore only plot values which are lighter than 
the strange-quark mass. The general picture here is sim- 
ilar to that seen for the degenerate-quark-mass pseu- 
doscalar and vector mesons: Again, perhaps partly due 
to somewhat large statistical errors for the excited states, 
the results for the pseudoscalar mesons appear consistent 
with the experimental values. The same also applies to 
the vector mesons on the fine lattices, while the coarse- 
lattice excited vector results appear systematically high. 



experimental value and, given the consistency of the two 
sets (within statistical errors), no large discretization ef- 
fects are apparent. 

In Fig. [51 we show the ground- and excited-state masses 
of the vector, pseudovector, and scalar mesons. For the 
vector mesons on the coarse lattices, only one excited 
state is extracted (another, slightly lighter effective mass 
plateau was seen, but since the eigenvectors were not sta- 
ble over the same region, no fits were performed). This 
state either corresponds to the p(1700) or it suffers ap- 
parently significant discretization effects. On the fine lat- 
tices, two such excited states are resolved which appear 
consistent with the physical states; the statistical errors, 
however, are larger, possibly hiding any residual system- 
atic effects. 

The pseudovector and scalar mesons present more dif- 
ficulties. Fewer effective mass plateaus can be found for 
excited states; for the b\ and scalar mesons, they are vir- 
tually non-existent. There are already significant prob- 
lems for the ground states. At small quark masses, it be- 
comes obvious that the pseudovectors suffer large quench- 
ing and/or finite- volume effects, causing an apparent en- 
hancement of the mass. 

The upward curvature seen in the ground-state vectors 
(at least on the coarse lattices) and pseudovectors may 
be partly explained by quenching effects. Quenched chi- 
ral perturbation theory for vector mesons |22j predicts a 
negative M„ term, which could certainly cause the ob- 
served enhancement at low quark masses. 

For the scalar meson plot similar curvature is seen, 



C. Comparison with results on smaller lattices 

Figure[7]displays two pseudoscalar effective mass plots, 
one from one of the current data sets (16 3 x 32) and one 
from a previous set (12 3 x 24) at the same lattice spacing 
(see Refs. |HE3)- O ne can see that previously chosen 
fit ranges (for the 12 3 lattices) may have been started at 
too small a value of t, apparently enhancing the resulting 
masses when compared to our current results. Since the 
physical spatial volume has been increased from 1.8 to 2.4 
fm (the latter perhaps still being too small for excited- 
state spectroscopy; see Ref. [25|), one may be tempted to 
simply put this down as a finite-volume effect. However, 
since the later effective masses on the smaller lattices are 
consistent with the plateau found on the bigger ones, it is 
not easy to say whether this state is displaying effects due 
to the finite volume. The enhancement may also be due 
to higher state corrections, which still contribute since 
the minimum timeslice considered in the fit is too small. 
With better statistics on the 12 3 lattices (to make up 
for the smaller number of spatial sites), we may end up 
finding a similar plateau to that seen on the 16 3 lattices. 
Similar effects are seen for the vector mesons. We suspect 
that this poorer level of statistics for the 12 3 lattices, and 
the subsequent choice of fit ranges, may account for most 
of the relative enhancement of our earlier results for the 
excited meson masses. 



6 



p(1700) C] jl f 









-1 1 


>■ 




1 % 


= s > 




n ' - 

Ms 






^10(1680) 





p(1450) 
p(770) 



,*0 * 

1020) 



- 0::H + )a=0.12fm->i i«-a=0.15fm(on) 





FIG. 5: Ground- and excited-state meson masses vs. M% for both lattice spacings. Vector (V), pseudovector (PV), and scalar 
(SC) mesons appear. The quark masses are degenerate. All quantities are in units of the Sommer parameter, ro- The diamonds 
represent the experimental points and the vertical lines mark the values of r%M^ corresponding to the physical strange quark 
mass. 



D. Chiral extrapolations of fine lattice results 

Although we lack a third lattice spacing which would 
allow us to try to work in the continuum, we neverthe- 
less attempt chiral extrapolations using the results from 
just our fine lattices (see Fig. |SJ|. For the most part, 
we use simple linear extrapolations (interpolations for <j>) 
in TqM% and, when doing so, avoid the low quark mass 
region (tqM% < 2.5) for the pseudovectors since these 
suffer obviously large systematics there. 

Assuming quenched chiral effects similar to those ex- 
perienced by vector mesons |22j. we also allow for a r^M^ 



term (obviously negative) in the fits for the pseudovec- 
tor extrapolations and include all points. We also try a 
quadratic fit in TqM% for the scalars. These additional 
extrapolations are plotted as the bursts on the right of 
the column for each of these states. The right hand side 
of the plot shows the persistent difficulties of simulat- 
ing pseudovector mesons in the quenched approximation. 
Dynamical simulations of these states appear not to suf- 
fer the same systematic enhancements |26l |27j . The ex- 
trapolations of the scalar appear more consistent with 
the ao(1450), a finding which appears consistent with 
unquenched results of simple ud scalars [23, 0] • 
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FIG. 6: Ground- and excited-state meson masses vs. M% for both lattice spacings. Both pseudoscalar (left) and vector (right) 
mesons appear. One quark mass is fixed to the physical strange quark mass. All quantities are in units of the Sommer 
parameter, ro. The diamonds represent the experimental points and the vertical lines mark the values of r%M^ corresponding 
to the physical strange quark mass. 
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FIG. 7: Ground- and excited-state effective masses for the 
pseudoscalars. Results for the 12 3 x 24 (see Refs. 0, Hc|p 
and 16 3 x 32 coarse lattices (a ~ 0.15 fm) are shown, along 
with horizontal lines displaying M ± a M from fits to the cor- 
responding time intervals. 



To the left in Fig. [S] we see that the results for the 
first-excited pseudoscalar and vector states agree with 
the experimental values (and dynamical calculations of 
excited pseudoscalars [22]]), while the ground-state vec- 
tors (except for the 4>) exhibit residual systematic effects, 
a common problem in the quenched approximation. 



Given the fact that our results include (possibly large) 
systematic errors (quenching, finite volume, and chiral ef- 
fects) which are not presently well understood, we point 
out that the correspondence between our masses and the 
physical ones in Fig. [S] should be viewed with some cau- 
tion. Nevertheless, despite these possible pitfalls, it is 
remarkable that the mass splittings appear to be of the 
right size. 



IV. SUMMARY 

We have presented ground- and excited-state meson 
masses from quenched lattice calculations using chirally 
improved quarks. Using a collection of linearly indepen- 
dent interpolating fields and the correlator matrix tech- 
nique, we have been able to clearly separate the different 
mass eigenstates. Although it must be admitted that our 
statistical errors are still quite large, the masses we find 
for excited pseudoscalars are consistent with experimen- 
tal values. The same is true for the excited vector mesons 
on our fine lattices. Results for the vectors on the coarser 
lattices suggest either that we have "missed" the first 
excitation (slightly lighter effective mass plateaus were 
seen in another eigenvalue, but without corresponding 
plateaus for the eigenvector) or that this state experi- 
ences significant discretization effects. The pseudovector 
and scalar mesons proved more difficult to handle. We 
could obtain few results for the excited states, virtually 
none below the strange quark mass. At low quark masses, 
already the ground states have problems, displaying sig- 
nificant quenching and/or finite-volume effects. Ulti- 
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mately, these problems will only be overcome via simula- 
tions with better interpolators, larger physical volumes, 
and dynamical fermions. 
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FIG. 8: Chirally extrapolated (interpolated for cj>) results 
from our fine lattices (20 3 x 32, a — 0.119 fm). Most extrap- 
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